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universe, which may be occurred in the early universe or late time, is studied. 



^ i' By considering the logarithmic correction to the energy density, we study the behavior 

C ' of Hubble parameter in the holographic dark energy model. We assume that the universe 
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I. INTRODUCTION 

To explain the cosmic accelerated expansion of the universe 11, and motivated by the holo- 
graphic principle [4], a model of dark energy has been proposed [SHOl] which has been tested and 
constrained by various astronomical observations |10|. This proposal is generically known as the 
"Holographic Dark Energy" (HDE). Its definition is originally extracted from the entropy-area re- 
lation which depends on the theory of gravity. In the thermodynamics of black hole horizons, there 
is a maximum entropy in a box of length L, commonly termed, the Bekenstein-Hawking entropy 
bound, S ~ rripL^, which scales as the area of the box j4 ~ L^ rather than the volume V ^ L^. 
Here rrip = (SvrG)"^ is the reduced Planck mass. In this context, Cohen et al. |3| proposed that 
in quantum field theory a short distance cutoff A is related to a long distance cutoff L due to the 
limit set by formation of a black hole, which results in an upper bound on the zero-point energy 
density. In line with this suggestion, Hsu and Li [J, |5|] argued that this energy density could be 
viewed as the holographic dark energy density satisfying pd = 'inP'm'i/L'^, where L is the size of a 
region which provides an IR cut-off, and the numerical constant 3n^ is introduced for convenience. 

It is essential to notice that in the literature, various scenarios of HDE have been studied via 
considering different system's IR cutoff. In the absence of interaction between dark matter and 



, Lila 



dark energy in flat universe, Li (5[ discussed three choices for the length scale L which is supposed 
to provide an IR cutoff. The first choice is the Hubble radius, L = H~^ J], which leads to a 
wrong equation of state for dark energy (w = 0), namely that for dust. The second option is the 
particle horizon radius. In this case it is impossible to obtain an accelerated expansion. Only the 
third choice, the identification of L with the radius of the future event horizon gives the desired 
result, namely a sufficiently negative equation of state to obtain an accelerated universe. However, 
as soon as an interaction between dark energy and dark matter is taken into account, the first 
choice, L = H~^, in fiat universe, can simultaneously drive accelerated expansion and solve the 
coincidence problem 111. It was also demonstrated that in the presence of an interaction, in a 
non-fiat universe, the natural choice for IR cutoff could be the apparent horizon radius |12l |. 

As we earlier mentioned that the black hole entropy S plays a central role in the derivation 
of HDE. Indeed, the definition and derivation of holographic energy density {pd = "in^mi/ L?') 
depends on the entropy-area relationship S ~ j4 ~ L^ of black holes in Einsteins gravity, where 
A '-^ L? represents the area of the horizon. However, this definition can be modified from the 
inclusion of the effects of thermal fluctuations around equilibrium |l5| | , quantum fluctuations [16|] , 



or by considering the loop quantum gravity (LQG) 



17|, all leading almost to the same result. 



The corrections provided to the entropy-area relationship leads to the curvature correction in the 



Einstein-Hilbert action and vice versa 



13(. The corrected entropy takes the form 



a 



where 7 and 5 are diniensionless constants of order unity. The exact values of these constants 
are not yet determined and still debatable. Taking the corrected entropy-area relation ([T|) into 



account, the energy density of the HDE will be modified as well. On this basis, Wei [1^ proposed 
the energy density of the so-called "entropy-corrected holographic dark energy" (ECHDE) in the 
form 

Pd = Sn^mlL"^ + 7L-^ ln{mlL'^) + 6L~^. (2) 

In the special case 7 = (5 = 0, the above equation yields the well-known holographic energy density. 

To solve some essential problems in standard cosmology, it is believed that there was also 

another stage of accelerated expansion in the early universe known as inflation [191 ]. In the same 

way that some models, proposed to explain the inflation( such as a scalar field in slow roll model 



2011). have been used to explain the present acceleration of the universe, one can examine different 
models of dark energy, to explain the evolution of the universe in the inflation era. 

In this paper we study all possible behaviors of the Hubble parameter in a universe dominated 
by (ECHDE) and a barotopic matter and study some consequences of this model in accelerated 
expansion in late time and also in the inflation era. 

The plan of the paper is as follows: In section 2 we construct a cosmological model and derive 
some useful expressions for our further uses. In section 3, we obtain different possible behaviors 
of Hubble parameter in the interacting ECHDE model. In section 4 we discuss the possible im- 
plications of this model in the accelerated expansion of the universe in the early stage and in 
the present epoch separately(we do not intend to study the whole history of the universe in a 
same framework). We also check the conditions under which the universe will undergo multiple 
acceleration-deceleration phases. Finally we conclude this paper in the last section. We use units 
h = G = c = kB = l- 



II. THE MODEL 

We consider a spatially flat Friedman Robertson walker (FRW) universe dominated by dark 
energy and matter (this can be dark matter, radiation and so on). The Friedman equations read 

H = -ATT{Pd + Pd + Pm), (3) 

where pd and pm are densities of dark energy and matter respectively. P^ is the pressure of dark 
energy and the Hubble parameter, H, is assumed to be differentiable. As we have noticed in the 
introduction, the corrections ([2]), can be obtained from the loop quantum gravity, as well as by 
the inclusion of the effects of thermal fluctuations around equilibrium , quantum fluctuations, or 
by considering charge or mass fluctuations. So it is safe to study this problem in the framework of 
FRW cosmology, using ordinary Friedman equations as was mentioned in [18l |. 

The matter and dark energy are allowed to exchange energy 21[ via the source term Q: 



Pd + 3H{wd + l)pd = -Q 

Pm + 3i7(l + Wm)Pm = Q- (4) 

Because of this interaction term, we have not the conservation of partial stress-energy tensors of 
matter and dark energy '■T^'matterY ~ ~'^(dark)- ^ ^' ^^ ^^ assumed that the matter component 
and dark energy have the same velocity which is the velocity of the whole fluid, V. We can write 

In the scalar fleld model of inflation, after the slow roll regime, the scalar fleld, whose energy 
density is pd, decays to radiation via a rapid coherent oscillation. The source term for this decay 
which allows the radiation creation from inflaton is taken as [22 1 



Q = aHpd, (6) 

where a is a constant. Also, in dark energy models different interactions between matter and dark 
energy are assumed. As the nature of dark energy has not yet been known, these interactions are 
taken from other models such as string theory and scalar tensor theory and so on, or as 

Q = ^H{pd + Pm) 

Q = rjHpm, (7) 



ly to alleviate the coincidence 



23| . To study the evolution of 



where (3 and r] are real constants, are proposed phenomenological 
problem and also to prevent the universe to undergo the big rip [z 
the universe we are obliged to make use of a specific interaction, which we choose a general form 

as y , 



Q = 3H{XmPm + ^dPd) (8) 

where Am and A^ are real constants. ([8]) reduces to ^ and d?]) for Am = 0, A^ = Am and A^ = 
respectively. Note that ([8]) is the same as the scalar 

1, 



Q = ^^.y.VZ, (AmTf^^,,,,) + A,T(7^,,) j , (9) 

written in the comoving frame. The equation of state (EoS) parameter of dark energy, tM^, is 
defined by Pd = WdPd and the (EoS) parameter of matter, < Wm = -f^, is assumed to be a 
constant, e.g. for cold dark matter we have Wm = and for radiation Wm = q- 

Different models have been proposed for dark energy, hereinafter we adopt the (ECHDE) model 



18 1 for which infrared cutoff is taken as L = ■^. In this model, the dark energy density may be 



expressed as 



"-^(l-^-f^^' + ^l- (^o) 



c, a and (3 are dimensionless real constants and their values are still in debates in the literature 



When j^ In(L^) + -^ ^ jd the model reduces to the ordinary holographic dark energy model. 
The correction terms are relevant in the early universe, and also in the late time provided that 
—2" becomes larger with respect to the present time. Note that in the units used in this paper the 

rrip 1—1 

reduced Planck mass is given by mp = (87r)''~^'^', so c^ is the same as v? in [l8|. The ratio of dark 
energy density to critical density is then 

na = c'^-aH^ln{H^) + PH^. (11) 

By construction we must have 

< c^ - aH"^ ln(i?^) + pH"^ < 1. (12) 

Time derivative of il^ is obtained as 

nd = -2HH~^{aH^ -Qd + c^)- (13) 



In terms of the (EoS) of the universe, 



T3l) can be written as 



From Q one can derive 



Pd + Pm 3 H^ 



nd = 3H{w + l)iaH^-nd + c^)- (15) 



r = 3Hr [wd-Wm + (Ad + rXm) ( ) ) , (16) 



where r = ^^. Using r = ,-, '' , and 

-U; = Wdi^d + "Wmf^m., (17) 

we obtain 

1 ^d ArfOrf ~ 

For a = /3 = 0, we have 0^ = c^ and u; = w^- If w^m = (e.g. when the matter is cold dark 
matter) and in the absence of interaction we get w = Wd = 0, implying that pd is the same as dark 
matter. This was the motivation of taking another infrared cutoff for the model in [5|]. By taking 
the interaction into account, the a = /3 = model can describe an accelerating universe with a 
non-dynamical 0,d corresponding to a scaling solution. 

From (|15p and (jlSp we obtain 

(Am - Arf - Wm. + l)^d- aH'^ -C^ -Xm + Wm ,-, ^. 

" = aH^-2nd + c- + l ' ^''^ 

which results in 

j:j _ ^H {Xd — Am + l)^d + Am — 1 ,„„s 

2 aH^-2nd + c'^ + l ' ^ ^ 

where Am = Am — Wm,. (j20p can be written as an autonomous first order differential equation 

where 

f{H) = {Xrn-l) + {Xd-Xm + l){c'' + (3H^-aH^ln{H^)) 

g{H) = -c^ + l + {a-2f3)H'^ +2aH^ln{H^). (22) 



Note that equation ()2ip requires that H and ah of the higher order time derivatives of H be zero 
at the time when H = 0, i.e. H = can occur only asymptotically. As a result H can not cross 
H = (note that at H = 0, H = 0). So we may assume H(t) > 0, Vt in the domain of validity of 
our model. We also assume that H is differentiable, therefore g{H) ^ and the sign of g{H) does 
not change. 

III. CLASSIFICATION OF THE HUBBLE PARAMETER BEHAVIORS 



Obtaining an analytical general solution to (j21|) is not possible. In this part, instead of solving 
(j2ip . by using some mathematical methods based on the properties of Lambert functions, we 
discuss and classify different possible behaviors of the Hubble parameter dictated by this model, 
in terms of its parameters (Am, A^, c^, a,/3). As we do not fix the parameters, various behaviors 
for the model are derived. Note that the results obtained in this section are general and do not 
necessarily match with our present or early universe. We will discuss this issue in the next section 
where as we will see only some cases in the classification have necessary (although not sufficient) 
conditions to describe the early and late time acceleration of our universe. 

For H + H^ > (<)0 the universe is accelerating (decelerating). For H > the universe is in 
super- acceleration (phantom) phase and for H < the accelerated universe is in the quintessence 
phase. Following the discussion in the last part of the previous section, the sign of H is unchange- 
able, hence the system is still in quintessence or phantom phase and quintessence-phantom crossing 
does not occur. 

By defining u = H^, we have 

H + H^ = ^u, (23) 

where F{u) := A + Bu + Cnln(n), and 

^ = 2(1 - c^) + 3p 

B = 2{a-2/3) + 3q 

C = 4a + 3s. (24) 

p = Am - 1 + (Ad - Am + l)c^ 

q = (Arf-Am + l)/3 

s = -(Arf-Am + l)a. (25) 



We have defined 



So / in ()22p can be written as f{u) = p + qu + su In(ti). 

To obtain the number of critical (fixed) points of the equation (I2ip . and to get some insights 
about the behavior of the system, we must find the number of zeroes and the behavior of f{H) 
in terms of the parameters of the model. In the same manner, we must also study the behavior 
of F[u). f and g in ()22|) and F in ()23p have the functional form K{u) := a + bu + culn{u). The 
general behavior of K in terms of its parameters is discussed in detail in the appendix. 

G{H) (or H) in ()2ip has a zero at ff = and at most two other zeroes at Hi = -y/ui and 



25( 1 . as follows (the real branches 



H2 = \/u2 determined in terms of Lambert W function: W{x) 
of W{x) are denoted by Wq and W-i): 

ui = exp ( W-i ( — exp ( - ) ) 

U2 = exp {W, (-^ exp (|)) - ^) . (26) 

These are fixed points of the equation (j2ip . hence they can not be crossed. We consider three 
domains: Di = (0,iJi), D2 = (i?i,i^2) and D^ = (i^2,oo). If 3t, such that H{t) £ Di then H{t) 
will restricted to Di. The sign of H depends on the sign of G in Di, if G{H) < (>)0 for H € Di, 
then H < (>)0. The stability of the model at critical points depends also on the sign oi g{H) in the 
region where H is restricted. We assume that H is differentiable (note that by construction g{H) 
is well defined: \/H € ?R.,g{H) € K). Hence g{H) has no roots in the region where H is restricted, 
and therefore its sign does not change. For simplicity we only study the cases with g{H) > 0. The 
cases with g{H) < can be treated in the same manner. Indeed, ^{Hi) < implies that H > 0. 
This requires that -^{Hi) and g{Hi) have different signs {Hi is the critical point). For a g{Hi) 
with opposite sign, H < 0. Note that g{H) > is the only physical choice in regions including 
H = 0: limi^^o 9{H) = 1 - 0^(77 = 0) > 0. 

Using the results obtained in the appendix and the above arguments, we can classify behaviors 
of the Hubble parameter as follows: 

For the very special case s = 0, a 7^ implies A,^ — A^, + 1 = 0. In this situation the only 
fixed point of the autonomous differential equation (j2ip is H = 0. If Am > 1, then H > which 
gives rise to a super accelerated expansion. For Xm < 1, H = becomes stable fixed point and for 
Am = 1 we obtain a de Sitter space time H = 0. 

According to the appendix we have the following possibilities (in situations where Hi = H2 we 
denote the root of / by i^i): 

(1,1) (s > 0,p > 0,(7 > —s (ln(|) + 1)): We have H{t) > 0, so i7 = is an unstable critical point 
and limt_^oo H = 00. 



(1,2) (s >0,p>0,q = -s (ln(f) + l)): If H{t) G iO,Hi), then H{t) > and limt^ocH{t) = 
Hi. If H{t) > Hi, H{t) > and limi_,oo H = oo. 

So (1,1) and (1,2) indicate that the expansion of the universe is super-accelerated. 

(1,3) (s >0,p>0,q< -s (hi(2) + l)): If H{t) G (0,iJi), then H{t) > corresponding the 
super-acceleration. For H{t) G {Hi,H2), we have H{t) < 0, which does not necessitate the 
acceleration of the expansion. For H{t) G {H2,oo), H{t) > 0. Here Hi is a stable critical point. 

(1,4) (s > 0,p < 0): If H{t) G (0,Fi), then H{t) < 0. For H{t) G (i?i,oo) we have H{t) > 0. 

The same analysis can be done for (II) situations (see the appendix): 

(II,l)(s < 0,p < 0,g < -s (ln(f) + l)): We have H{t) < so \imt^ooH{t) = 0. 

(11,2) (s <0,p<0,q = -s (ln(f) + l)): U H{t) G (0,iJi), then H{t) < and \iint^ooH{t) = 0. 
If H{t) > Hi, then H{t) < and limt_,oo H{t) = Hi. 

(11.3) {s<0,p<0,q> -s(ln(f) + l)): liH{t) G {0, Hi), then H{t) < 0. For ii'(t) G {Hi,H2), 
we have i:r(t) > and finally if H{t) G {H2,oo), then -ff(t) < 0. In this situation H2 is a stable 
critical point. 

(11.4) (s < 0,p > 0): If H{t) G (0,ii'i) then H{t) > 0. But if H{t) G (ii'1,00) then H{t) < 0. 
Here -ffi is a stable critical point. 

In the above, situations corresponding to -ff > are related to the phantom phase. We remind 
that in an open set {0,Hi), {Hi,H2) or {H2,oo), the sign of H does not change and, e.g. if the 
universe is in phantom phase {H > 0) in some time of inflation, it will remain in this phase as long 
as the model is valid. The same is true for H < 0. Note that the cases corresponding to H < 
do not necessitate an accelerated expansion. In this case to see whether there is an acceleration 
phase we must study H + H^. 

So let us examine H + H^. This expression, besides at H = 0, has at most two positive zeroes 
at H3 and H4^: 



2 , .., , A fB\\ B 



^3 ^exp(„,(-ie.p (§))-§). (27) 

As these are not fixed points, they can be crossed allowing consecutive acceleration deceleration 
phases, but note that these transitions may only occur when H < 0. 

Again we can classify generally the model as follows (to avoid any confusion, we use /// and 
IV instead of / and // used above, but note that [III,i), and {IV, i), again, correspond to {I,i) 
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and {II, i) in the appendix respectively, and in situations where H-^ = H/^ we denote the root by 

(111,1) and (III,2)(C>0,A>0,B > -C (hi(^) + l)): We have H + H'^{t) > 0. 

(III,3)(C>0,A>0,S < -C(ln(^) + l)): \i H{t) G {Q, H ^), ilien H + H'^{t) > 0. Yoi H{t) e 
{H^, Hi), we have H + H^{t) < 0. For H{t) G (^^4, 00), H + H^{t) > holds. 

(III,4)(C > 0, A < 0): If H{t) G (0, H3), then H + H^ < and for /7(i) > i^g: H + H^ > 0. 

The same analysis can be done for (IV) situations (corresponding to (II) in the appendix) 

(IV,1) and (IV,2) (C7 < 0,^ < 0, -B < -C7 (in(^) + l)): We hme H + H^ < 0. 

{W,3){C <0,A<0,B> -C7(ln(^) + 1)): li H{t) G (0,^/3), then H + H^ < 0. For H{t) G 
{H3, Hi), we have H + H'^ > and finally if H{t) G (i74, 00), then H + H'^ < 0. 

(IV,4) (C < 0,A > 0): If H{t) G (0,^/3) then H + H^ > 0. But if i/(t) G (i/3,00) then 
i/ + if 2 < 0. 

In summary by studying all various possible behaviors of the Hubble parameter as the solution 
of autonomous differential equation (|2ip . in (I-II, 1-4) the conditions leading to H > and H < 
and all possible late time behaviors of H are specified. In this model the transition from phantom 
to quintessence and vice versa are forbidden. This is due to the fact that the points where H = 0, 
are critical points of the theory. Although the conditions to have a phantom phase H > can 
be read from (1-11,1-4) but they don't elucidate the conditions required to have the quintessence 
phase. So by studying all possible behaviors of H + H^, conditions needed for accelerated and 
decelerated expansion were specified. 

IV. PHYSICAL DISCUSSION AND RESULTS 

In this part we try to study physical implications of our model in the late time and separately 
in the inflation era of our universe. 

A. Late time acceleration 

From I and II, we flnd out that all the possible late time solutions in this model are de Sitter 
space-time and the asymptotic values of the Hubble parameter may be the critical points: H = 0, 
H = Hi or H = H2 specified by (|26p . We remind that the system is restricted to the regions 
Di = (Hi, Hj) bounded by critical points and the sign of H is unchangeable in Z?j. So if i7 < (>)0 
in Di, then liuit^oo H = Hi{Hj). The model is stable at H^ where f{Hi) = provided that 
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q(h ) dJ{(-^i) ^ *-*' ^^^ ^ necessary condition for the system to tend to i/ = is ^L| < 0. For 
liiTii^oo H = 0,we have ^(0) = and ^ = ^ >'rn~i+{^d^~>^r.+i)c^ ^ q j^ ^^le absence of interaction 
this inequaUty is satisfied but the interaction may prevent the model to go asymptotically to H = 0, 
and instead, forces it to tend to Hi. 

In the present era we have —3- <^ 1, therefore -^ln(L^) + "^ ^ T2 cind the correction terms 
have not important role. So at first sight, it seems that the corrections are not relevant in the late 
time acceleration, but this is not true. To emphasize this via a simple example, let us take the 
ordinary holographic dark energy model. In the absence (or smallness) of the correction terms, 
the only critical point of (I2ip . is -ff = 0. If we ignore the interaction, we have H = —^H^ < and 
lim(_!,oo H{t) = and the universe tends to a static space time at late time, besides, H + H^ < 
and the acceleration does not occur. During this evolution Q^ has the constant value Jl^^ = c^. 
Here the correction terms will be also irrelevant in the future evolution of the universe. Now let us 
take the interaction into account and let ( (Am — 1) + jz^ ) < (this corresponds to p > cases 
in (I) and (II). Then H > and the correction terms will become relevant at late time. 

As in our epoch —3- <^ 1, the present time belongs to the region (0,i7i) and ii H < 0, then 
there is no need to correction terms to study the future evolution of the universe and we have 



li. If, 



limj_!.oo H{t) = 0, so the evolution of the universe can be explained in the same manner as 
instead, H > which is related to (1,1), (1,2), (1,3) or (11,4) the correction terms become important 
and the universe tends to a de Sitter space time characterized by : limi_j.oo H{t) = Hi. (1,1) implies 
limj_j.oo H = 00 and is excluded by[T2j So among various cases in our classification in the previous 
section we are left only with s > 0, p > 0, q < —s (ln(^ + 1)) and the region (0, Hi). 
The value of il^ at late time is 

lim na = , ^"Z" . (28) 

t^-oo Ad — Am + 1 

So the interaction, besides preventing the model to tend to H = 0, via the energy exchange, 
alleviates the coincidence problem. In this case it is worth noting that if a = /3 = 0, then 
lim(_j.oo H(t) = 00. So the correction terms prevent the Hubble parameter to become very large 
asymptotically. 

B. Inflation 

Among various possible acceleration phases, it seems that the phantom phases, reported through 
the situations I and II, are not consistent with a non eternal inflationary phase. Indeed as long 
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as the universe is dominated by (ECHDE) this inflationary phase continues, i.e. the system is 
restricted to the domain specified by critical points. So if the phantom phase occurs, can not 
be ceased. Although crossing the critical points, Hi and H2, is not possible but the system can 
cross ^^3 and H4 and therefore transition from acceleration to deceleration and its inverse are 
possible. To see this, as an example, let us take the case (111,3) where for H{t) G {H^, H4), we have 
H + H^ < therefore H < and H is decreasing. But H^ is not a fixed point and the Hubble 
parameter crosses H^ and enters in accelerating domain. 

To describe the temporary inflationary phase, we must consider H < situations in I and II 
and then investigate the cases reported in III and IV. If we expect that the inflation be ended, we 
must select the situations allowing the transition from acceleration to declaration phase. This is 
related to the presence of the source term Q, allowing the energy exchange between dark energy 
and matter. This is only allowed in (111,3) with Hmf. > H^, (III>4) with H-mf. > H-:^, and (IV, 3) 
with Hinf, £ {H^, H4), provided that the H belongs to cases in I and II where H < 0. Hi^f denotes 
the value of the Hubble parameter during inflation. At the end of inflation, H is determined by H^ 
(cases : (IV,3)and (III,4))or H4 (case (111,3)) determined by (f27|l . The corresponding values of fi^ 
can then be read from (jlip . In the above cases it is straightforward to see that (IV) corresponds 
to an inflation which is not past eternal. 

Until now, for the sake of generality we did not fix the values of parameters {a, /?, c^. Am, A^}. 
To be more specific, and as an illustration of our results, let us take the case (11,1) and (111,4) 
characterized by 

(Am - Arf - l)a < 

(Am-l)(l-c2) + Arfc2<0 (29) 

and 

(1 - 3Arf + 3Am)a > 

(3Arf - 3Am + l)c2 + 3Am - 1 < 0, (30) 

respectively. This situation describes a deceleration followed by an acceleration phase as can be seen 
from the diagram of H + H^ depicted in fig.([II for the optional choice {a = 4.18, /3 = —10.71, (P = 
0.7, Am, = 0, Arf = 0.1,} which satisfies ([2^ and ([5U|) (note that the same behavior is expected for 
all models whose parameters satisfy (j29p and (j30p ). The inflation ends at H^ = 0.02(~ O.Gmj,) 
and ^ restricts H'^ to H"^ < 0.18(0. 95m|>). 
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FIG. 1: H + H'^ smsi function oi u = H^, for a = 4.18, /3 = -10.71, A^ = 0.1, A™ = 0, c^ = 0.7 



It is also interesting to see if the model permits to have consecutive acceleration and deceler- 
ation phases. We remind that H + H^ has at most two zeroes besides H = 0. So in general, 
besides acceleration-deceleration or deceleration-acceleration transitions, it is also possible to have 
successive deceleration-acceleration phases in this model. Acceleration-deceleration-acceleration 
corresponds to (111,3) and deceleration-acceleration-deceleration corresponds to (IV, 3), provided 
that the zeroes oi H + H^ belongs to domain Di (specified by critical points) where the Hub- 
ble parameter is restricted. Now to illustrate this result let us take the case (11,1) and (111,3) 
characterized by (j29j) and : 



a(3A™ - 3Ad + 1) > 

3A^ - 1 + (3(Ad - A^) + 1) c^ > 

3(A^ - Ad + l)(a - /3) + 2/3 < a(3Ad - 3A, 
^^(3A„-l)(l-c2) + 3Arfc2 



l)x 



(31) 



a(3Am - 3Ad + 1) 

H + H^ is depicted in terms of n = H^, for the choice {a = 0.5, /3 = —0.5, A^ = Am = 0.103, c^ = 
0.7}, which satisfy (j29p and (j3ip . in fig. ([2]) (again note that the same qualitative behavior is expected 
for all models whose parameters satisfy (j29p and (j3ip ). The model has an acceleration phase for 
H"^ > 0.0218 (or m h = c = 1 units H"^ > 0.548m|,) a deceleration phase for 0.0150(= 0.377m^) < 
H^ < 0.0218(= 0.548m|>) and again an acceleration phase for < H^ < 0.0150(= 0.377m^). Note 
that (fT2]l restricts the model to H^ < 1.19 = (29.9mp). Models with double inflation (i.e. two 
stages of inflation) were studied and reported before in the literature |26l | . 

The scope of our discussion was only restricted to the study of possible accelerations of the 
universe in early and late times. However a realistic general model of inflation must also be 



capable to describe physical problems such as the reheating and the cosmological perturbations 



14 




0.005 0.01 COlV^ 0.02/ 0.025 0.03 



FIG. 2: iJ + i/2 as a function of u = H'^, for a = 0.5, /3 = -0.5, A^ = A™ = 0.103, c^ = 0.7 



and so on. 

In the inflationary era, the universe is dominated by p^ for which the pressure is negative and by 
giving rise to H + H^ > drives the inflation. In different models, there may be different scenarios 
for the reheating. As an example, in the scalar field model and in the slow roll approximation, 
after the infiation the scalar field decays to other particles during a rapid coherent oscillation in 



the bottom of the potential slope 22] • In our model and in the absence of scalar fields, there can 
be another possibilities for reheating: 

By taking Wm = 1/3 and Xm = 0, the model acts similar to warm infiation models where 
the infiaton decays to radiation (ultra-relativistic particles) during infiation [27|. This is due 
to the source term, SHX^Pd, which allows energy exchange between these components. In the 
infiationary era, the universe is dominated by pd, and this source term is significant in (jl|). In 
contrast, after the infiation, the contribution of pm becomes more significant and the universe 
becomes radiation dominated for H + 2H^ ~ which occurs in deceleration era. In this era the 
radiation creation from p^ is not significant. Note that the radiation dominated era occurs for 
H + 2H^ ~ 0. The corresponding value of the Hubble parameter can be obtained by solving the 
equation Ar + BrU+ Cruln{u) = 0, where Ar = 4(1 — c^) + 3p; B^ = 4(a — 2/3) + 3g and Cr = 8a + 3s. 

One can also assume another possibility for reheating by taking Wm ^ 4- This assumption 
may be valid provided that pm. decays to ultra-relativistic particles at the end of infiation, giving 
rise to the preheating or the reheating. However the study of this era requires that we consider 
the contribution of baryonic matter in Q and add the corresponding interaction terms to our 
equations. 

At the end it is worth to note that the above conditions posed on the parameters of the model, 
(q,/3. Ad, Am), in (III, 3), (111,4), and (IV,3) although are necessary conditions for transient infiation 
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but are not sufficient. For example for a given model with specified parameters the e-folds number 
A^ := In ( ^]^[p\ ) must be calculated. It is given by 

where as before u := H^, i and eiid denote the beginning and the end of the inflation and G{u) is 
given by (f2T]l . In the case (111,4); (111,3); and (IV,3) we have Uend = -fff; Uend = H^', and Uend = H^ 
respectively. In the case (IV, 3), H^ < Ui < H^. Note that in other two cases Uj € (uendi^) where 
V is the maximum value of u satisfying (J12p . For a viable model we must have 60 < A^. This put 
more constraints on the parameter of the model. However for a general model, analytically solving 
([32]) or obtaining a lower bound for it, are very complicated tasks. 



V. CONCLUSION 

In this paper we considered a spatially flat FRW universe dominated by (ECHDE) and a 
barotopic matter interacting via a more general source term with respect to other papers in this 
subject. We took the apparent horizon as the infrared cutoff (this is the more natural choice 
adopted in the literature). Considering the effects of thermal fluctuations around equilibrium, 
quantum fluctuations, or charge and mass fluctuations, modifies the entropy attributed to the 
horizon. We applied these corrections to the apparent horizon entropy and achieved to obtain an 
autonomous differential equation for the Hubble parameter. Then we obtained the critical points 
of the model and classified the behavior of the system in terms of the parameters of the interaction 
and (ECDHE) . For this purpose we used algebraic features of the autonomous differential equation 
and LambertW functions. 

Although in the present time the corrections are marginal but they may play an important 
role in the early and late times. We deduced that the correction terms may force the universe to 
tend to a de Sitter space-time at late time. We obtained the possible ultimate value of the Hubble 
parameter and also derived the corresponding dark energy density. We showed that the coincidence 
problem is alleviated in this model. 

In addition we studied some necessary (although not sufficient) conditions for the model to 
describe the acceleration phase in inflation era. The inflation was assumed to be transient and the 
possible values of the Hubble parameter at the end of inflation were derived. 

Appendix 
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In this part we study the behavior and the properties of the roots of K(u) := a + bu + cu ln{u). 
y^' has and only has a root at u = exp (— ^-'-^), so following the Rolle's theorem K{u) has at 
most two positive roots which we denote ui and U2 > ui. 

We have also -7-3- (S) = §• Generally our model can be classified into three classes: c = 0, I: 
c > and II: c < 0. 

In the case I: 

1) for -ftr(O) > and K{u) > 0, K has no roots and K{u) is always positive. 

2) For K{{)) > and K{u) = 0, -ftT has only one root and K[u) is positive. 

3) For ivr(O) > 0, and K{u) < 0, K has two roots: < ui < U2, and K{u > U2) > 0, 
K{ui < u < M2) < and K{0 < u < m) > 0. In this case ^{ui) < and ^{m) > 0. 

4) For K{0) < and K{u) < 0, K has only one non zero root ui. We have also K{u > ui) > 0, 
K{0 <u<ui) <0 and ^(ni) > 0. 

The cases (1,1), (1,2), (1,3), and (1,4) correspond to 

c > 0, a > 0, 6 > -c Mn - + 1 

c > 0, a > 0, 6 = -c (in - + 1 

c > 0, a > 0, 6 < — c f In — I" 1 ) > and; 

0, a < 0, (33) 

respectively. 
In the case II: 

1) for K{0) < and K{u) < 0, K has no roots and K{u) is always negative. 

2) For K(0) < and K{u) = 0, K has only one root and K{u) is always negative. 

3) For K{0) < 0, and K{u) > 0, K has two roots: < ui < U2, and K{u > U2) < 0, 
K(ui < u < M2) > and K{0 < u < m) < 0. In this case ^{ui) > and ^{m) < 0. 

4) For K{0) > and i^(n) > 0, K may have only one nonzero root ui, and K{u > ui) < 0, 
ir(0 < u < ni) > and ^(ui) < 0. 

The cases (11,1), (11,2), (11,3), and (11,4) correspond to 

c < 0, a < 0, 6 < -c Mn - + 1 

c < 0, a < 0, 6 = -c Mn - + 1 j ; 

c < 0, a < 0, 6 > — c f In — I" 1 ) ) o^"^; 

c < 0, a > (34) 
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respectively. 

The numbers of the roots of K can also be explained in terms of Lambert W function: W{x) 



25l |. in a more straightforward way. The real branches of W{x) are denoted by Wq and W-i. For 
real x, if — -^ < x < there are two possible real values for W{x): —1 < Wq{x), and W-\{x) < —1. 
We have also Wo{-^) = VF_i(-i) = -1. 

For ^ec < (the cases 4 in I and II), the solution of the equation K = in terms of Lambert 



W function is 



u = exp(l^o(-^exp(^]]-^) (35) 



For - < -ec, K has no real roots (the cases 1 in I and II), and for < -e^ < -, K has two roots 
(the cases 3 in I and II): 

ui = exp I W-i I — exp I - I I 

„, = exp(M/„(-^exp(^))-^). (36) 

For - = -e~, ui = U2 corresponding to the cases 2 in I and II. 
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